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Abstract 

The swimming of a pair of spherical bladders that change their 
volumes and mutual distance is efficient at low Reynolds numbers and 
is superior to other models of artificial swimmers. The change of shape 
resembles the wriggling motion known as metaboly of certain protozoa. 

Swimming at low Reynolds numbers can be remote from common intu- 
ition because of the absence of inertia pp. In fact, even the direction of 
swimming may be hard to foretell j2J. At the same time, and not unrelated 
to this, it does not require elaborate designs: Any stroke that is not self- 
retracing will, generically, lead to some swimming [3]. A simple model that 
illustrates these features is the three linked spheres |3|, Fig. [T] (right), that 
swim by manipulating the distances between neighboring spheres. The 
swimming stroke is a closed, area enclosing, path in the l\ — i<i plane. An- 
other mechanical model that has actually been built is Purcell's two hinge 
model 

Swimming efficiently is an issue for artificial micro-swimmers [5]. As 
we have been cautioned by Purcell not to trust common intuition at low 
Reynolds numbers [2], one may worry that efficient swimming may involve 
unusual and nonintuitive swimming styles. The aim of this letter is to give an 
example of an elementary and fairly intuitive swimmer that is also remarkably 
efficient provided it is allowed to make large strokes. 

The swimmer is made of two spherical bladders, Fig. ^ (left). The blad- 
ders are elastic bodies which impose no-slip boundary conditions. The device 
swims by cyclically changing the distance between the bladders and their rel- 
ative volumes. For the sake of simplicity and concreteness we assume that 
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their total volume, Vq, is conserved. The swimming stroke is a closed path in 
the v — £ plane where v is the volume of, say, the left sphere and £ the distance 
between them. We shall make the further simplifying assumption that the 
viscosity of the fluid contained in the bladders is negligible compared with 
the viscosity of the ambient fluid. For reasons that shall become clear below 
we call the swimmer pushmepullyou. 

Like the three linked spheres, pushmepullyou is mathematically elemen- 
tary only in the limit that the distance between the spheres is large, i.e. 
when S{ = di/l <C 1. (eij stands for the radii of the two spheres and £ for 
the distances between the spheres.) We assume that the Reynolds num- 
ber R = pav/fi <C 1, and that the distance £ is not too large: £v <C fi/p. 
The second assumption is not essential and is made for simplicity only. (To 
treat large £ one needs to replace the Stokes solution, Eq. (JTTJ), by the more 
complicated, but still elementary, Oseen-Lamb solution 

Pushmepullyou is simpler than the three linked spheres: It involves two 
spheres rather than three; it is more intuitive and is easier to solve mathe- 
matically. It also swims a larger distance per stroke and is considerably more 
efficient |S] . If large strokes are allowed, it can even outperform conventional 
models of biological swimmers that swim by beating a flagellum j^j. If only 
small strokes are allowed then pushmepullyou, like all squirmers j^j, becomes 
rather inefficient. 

The swimming velocity is defined by X = (U\ + U^jl where Ui are the 
velocities of the centers of the two spheres. To solve a swimming problem 
one needs to find the (linear) relation between the (differential) displacement 
dX, and the (differential) controls (d£,dv). This relation, as we shall show, 
takes the form: 

2dX = ^1^1 d £ +-L dv, (1) 
ai + a 2 2n£ 2 ' v ; 

where a 1; a 2 are the radii of the left and right spheres respectively and v is the 
volume of the left bladder. dX stresses that the differential displacement does 
not integrate to a function X(£,v). Rather, the displacement X(j) depends 
on the stroke 7, defined as a closed path in £ — v plane. The first term says 
that increasing £ leads to swimming in the direction of the small sphere. It 
can be interpreted physically as the statement that the larger sphere acts as 
an anchor while the smaller sphere does most of the motion when the "piston" 
£ is extended. The second term says that when £ is held fixed, the swimming 
is in the direction of the contracting sphere: The expanding sphere acts as 
a source pushing away the shrinking sphere which acts as a sink to pull the 
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Figure 1: Mue snapshots of the pushmepullyou swimming stroke (left) and 
the corresponding strokes of the three linked spheres (right). Both figures are 
schematic. After a full cycle the swimmers resume their original shape but 
are displaced to the right. Pushmepullyou is both more intuitive and more 
efficient than the three linked spheres. 

expanding sphere. This is why the swimmer is dubbed pushmepullyou. 

To gain further insight consider the special case of small strokes near equal 
bi-spheres. Using Eq. (JTJ) one finds, dropping sub-leading terms in = Oijt. 

5X= -dlogvAdi (2) 
6 

The distance covered in one stroke scales like the area in \ogv — £ plane. Note 
that the swimming distance does not scale to zero with e, when the spheres 
are far apart. This is in contrast with the three linked spheres where the 
swimming distance of one stroke is proportional to e. For a small cycle in 
the i\ — £2 plane Najafi et. al. find for a symmetric swimmer (Eq. (11) in 

m 

5X = 0.7e d\og£ 2 A d£\ (3) 

When the swimmer is elementary, (=when e is small), it is also poor. 

Consider now a large stroke associated with the closed rectangular path 
enclosing the box £ s < £ < £l, v s < v±, v% < vl = fo — v s , where v 1 = v and 
t>2 are, respectively, the volumes of the left and right bladders. If a s C 
then from Eq. (JTJ), -^(7) is essentially £l — £ s - 

X(7)=(^0 (h-£ s )(l + 0(e 3 )) (4) 
\ol + a s J 
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This says that the distance covered in one stroke is of the order of the size 
of the swimmer, i.e. the distance between the balls t. 

Certain protozoa and species of Euglena perform a wriggling motion 
known as metaboly where, like pushmepullyou, body fluids are transferred 
from a large spheroid to a small spheroid pUj- Metaboly is, at present not 
well understood and while some suggest that it plays a role in feeding oth- 
ers argue that it is relevant to locomotion ^T]. The pushmepullyou model 
shows that at least as far as fluid dynamics is concerned, metaboly is a vi- 
able method of locomotion. Racing tests made by R. Triemer ^2] show 
that Euglenoids swim 1-1.5 their body length per stroke, in agreement with 
Eq. (@J) for reasonable choices of stroke parameters. Since Euglena resemble 
deformed pears — for which there is no known solution to the flow equa- 
tions — Pushmepullyou is, at best, a biological over-simplification. It has 
the virtue that it admits complete analysis. 

The second step in solving a swimming problem is to compute the power P 
needed to propel the swimmer. By general principles, P is a quadratic form 
in the velocities in the control space and is proportional to the (ambient) 
viscosity fi. The problem is to find this quadratic form explicitly. If the 
viscosity of the fluid inside the bladders is negligible, one finds that in order 
to drive the controls £ and v, Pushmepullyou needs to invest the power 

1 + iV , # + 2. (i + i-V (5) 



Note that the dissipation associated with £, is dictated by the small sphere 
and decreases as the radius of the small sphere shrinks. ( The radius can 
not get arbitrarily small and must remain much larger than the atomic scale 
for Stokes equations to hold.) The moral of this is that pushing the small 
sphere is frugal. The dissipation associated with v is also dictated by the 
small sphere. However, in this case, dilating a small sphere is expensive. 

The drag coefficient is a natural measure to compare different swimmers. 
It measures the energy dissipated in swimming a fixed distance at fixed speed. 
(One can always decrease the dissipation by swimming more slowly.) Let r 
denote the stroke period. The drag is formally defined by [§J E3] : 

r n Pdt 

X{pf) is the swimming distance of the stroke 7. The smaller 5 the more 
efficient the swimmer. 5 has the dimension of length (in three dimensions) 
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and is normalized so that dragging of a sphere of radius a with an external 
force has 5 = a. 

To compute the dissipation for the rectangular path we need to choose 
rates for traversing it. The optimal rates are constant on each leg provided 
the coordinates are chosen as (£, arcsin y^)- This can be seen from the 

fact that if we define x = arcsin ./—, then Avni 2 = { — + — ) v 2 and the 

Lagrangian associated with Eq. (JSJ is quadratic in (£, x) with constant coef- 
ficients, like the ordinary kinetic Lagrangian of non relativistic mechanics. It 
is a common fact that the optimal path of such a Lagrangian has constant 
speed. 

^From Eq. (JHJ) we find, provided also l\ 3> i 2 ., £lI&s 3> \JvlJv s 

Pdt^ 2 -^(l + o(e 2 ^^)), T e + T v = r/2 (7) 



67T/X J T t \ \ v s T t 

where Tt (T v ) is the time for traversing the horizontal (vertical) leg. (Here e 2 
is actually (a s /£L) 2 rather then the much larger (ai/is) 2 - Also note that the 
second term in Eq. (jSJ) contributed 0{vl/Tz) rather then 0{y 2 L / {y s T()) as one 
may have expected from Eq. (jSJ) which is dominated by the small volume.) 
The optimal strategy, in this range of parameters, is to spend most of the 
stroke's time on extending £. By Eqs. (I6I4I7|) this gives the drag 

5 w 4a s (8) 

where a s is the radius of the small bladder. This allows for the transport of a 
large sphere with the drag determined by the small sphere. To beat dragging, 
we need a s = a/4, which means that most of the volume, 63/64, must be 
shuttled between the two bladders in each stroke. 

It is instructive to compare Pushmepullyou with the swimming efficiency 
of models of (spherical) micro-organisms that swim by beating flagella. These 
have been extensively studied by the school of Lighthill and Taylor IT4*] 
where one finds 5 > 100 a. This is much worse than dragging. (We could 
not find estimates for the efficiency 5 for swimming by ciliary motion [To] , 
but we expect that they are rather poor, as for other squirmers jU].) For 
models of bacteria that swim by propagating longitudinal waves along their 
surfaces Stone and Samuel [12] established the (theoretical) lower bound 
5 > |a. (Actual models of squirmers do much worse than the bound.) If 
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the pushmepullyou swimmer is allowed to make large strokes, it can beat the 
efficiency of all of the above. 

Eqs. fl5!8|) do not strictly apply to metaboly because the viscosity of the 
fluid inside the organism can not be neglected and presumably dominates 
the dissipation. Euglena are not as efficient as Pushmepullyou. 

It is likely that some artificial micro- swimmers will be constrained to make 
only small (relative) strokes. Small strokes necessarily lead to large drag [0] , 
but it is still interesting to see how large. Suppose 5\og£ ~ <51ogt>, a\ ~ a 2 - 
The dissipation in one stroke is then 

^From Eq. (J2J) and noting that T e = ~t, one finds 

72 
(b logi;)^ 

We shall now outline how the key results, Eqs. (|1|5|) . are derived. The flow 
around a pair of spheres is a classical problem in fluid dynamics which has 
been extensively studied [13 HI] . We could have borrowed from the general 
results, e.g. in [TE], and adapt them to the case at hand. However, it is 
both simpler and more instructive to start from scratch: The classical Stokes 
solution \7[ describing the flow around a single sphere of radius a dragged by 
a force / and, in addition, dilated at rate v 

nu(x; a, /, v) = - J- ( ( 3 + Q f+(l- 3(/ ■ x)x] + ~ 2 x. (11) 
bfi\x\ \\ x A ) \ x z J J x A 

u(x; a, f, v) is the velocity field at a position x from the center of the sphere. 
The left term is the known Stokes solution. (A Stokeslet, [7], is defined as 
the Stokes solution for a = 0.) The term on the right is a source term. 

Since Stokes equations are linear, a superposition of the solutions for two 
dilating spheres is a solution of the differential equations. However, it does 
not quite satisfy the no-slip boundary condition on the two spheres: There is 
an error of order e. The superposition is therefore an approximate solution 
provided the two spheres are far apart. 

The (approximate) solution determine the velocities Ui of the centers of 
the two spheres: 

Ut = u(aj; a,, (-)V, 0) + u((-y£f; a v (-)*/, (-)V), i + j G {1, 2} (12) 
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The first term on the right describes how each sphere moves relative to the 
fluid according to Stokes law as a result of the force / acting on it. The 
second term (which is typically smaller) describes the velocity of the fluid 
surrounding the sphere (at distances 3> a but <C i) as a result of the move- 
ment of the other sphere. By symmetry, the net velocities of the two sphere 
and the net forces on them are parallel to the axis connecting the centers of 
the two spheres, and can be taken as scalars. To leading order in e Eq. (112)) 
reduces to 

^u i = {-yL(l--L\ + ^ (13) 



/ = -6tt// [- + -) I (14) 



fi \3ai 2£J 2£ 2 
Using £ = —U\ + U 2 gives the force in the rod 

y 

Dropping sub-leading terms in e gives Eq. (fTj). 

We now turn to Eq. Consider first the case v = 0. The power supplied 
by the rod is —f{Ui — U\) = —fi which gives the first term. Now consider 
the case i = 0. The stress on the surface of the expanding sphere is given by 



„=-5!*f^=4 (15) 
47T V x 2 J na 3 



The power requisite to expand one sphere is then 

4ira 2 ad = av = ^(v) 2 (16) 

Since there are two spheres, this give the second term in Eq. (pj. 

There are no mixed terms in the dissipation proportional to Iv. This can 
be seen from the following argument. To the leading order in e°, which is all 
we care about, the metric must be independent of i, (see Eq. Sending 
i —>■ — i is equivalent to exchanging the two spheres. This can not affect the 
dissipation and hence the metric must be even function of I. In particular, 
there can not be a term vl in the metric. This completes the proof of Eq. (j^J). 
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